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The active region of a ballistic nanostructure is an open quantum-mechanical system, whose 
nonunitary evolution (decoherence) towards a nonequilibrium steady state is determined by carrier 
injection from the contacts. The purpose of this paper is to provide a simple theoretical descrip- 
tion of the contact-induced decoherence in ballistic nanostructures, which is established within the 
framework of the open systems theory. The active region's evolution in the presence of contacts 
is generally non-Markovian. However, if the contacts' energy relaxation due to electron-electron 
scattering is sufficiently fast, then the contacts can be considered memoryless on timescales coars- 
ened over their energy relaxation time, and the evolution of the current-limiting active region can 
be considered Markovian. Therefore, we first derive a general Markovian map in the presence of 
a memoryless environment, by coarse-graining the exact short-time non-Markovian dynamics of an 
abstract open system over the environment memory-loss time, and we give the requirements for the 
validity of this map. We then introduce a model contact-active region interaction that describes 
carrier injection from the contacts for a generic two-terminal ballistic nanostructure. Starting from 
this model interaction and using the Markovian dynamics derived by coarse-graining over the effec- 
tive memory-loss time of the contacts, we derive the formulas for the nonequilibrium steady-state 
distribution functions of the forward and backward propagating states in the nanostructure's active 
region. On the example of a double-barrier tunneling structure, the present approach yields an 
I-V curve that shows all the prominent resonant features. We address the relationship between the 
present approach and the Landauer-Biittiker formalism, and also briefiy discuss the inclusion of 
scattering. 
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I. INTRODUCTION 

In a nanoscale, quasiballistic electronic structure un- 
der bias, relaxation towards a steady state cannot 
be described by the semiclassical Boltzmann trans- 
port equatioi*^, because the structure's active region 
is typically smaller than the carrier mean free path 
and efficient scattering no longer governs relaxation. 
Rather, the nanostructure's active region behaves as an 
open quantum-mechanical systemf^i^ exchanging parti- 
cles with the reservoirs of charge (usually referred to as 
leads or contacts). In the absence of scattering within 
the active region, the coupling of the active region to 
the contacts is the cause of its nonunitary evolution 
(decoherence) towards a nonequilibrium steady state, 
and the importance of this coupling has become well- 
recognized in quantum transport studies. The descrip- 
tion and manipulation of the contact-induced decoher- 
ence are presently of great importance not only in quan- 
tum transport studies ^n^n^nL^ ^^^i ^ ^ i but also in the the- 
ory of measurement^^ and quantum information^^ 

The purpose of this paper is to provide a simple 
description of the nonunitary evolution of a ballis- 
tic nanostructure's active region due to the injection 
of carriers from the contacts. Carrier injection from 
the contacts into the active region is traditionally de- 
scribed by either an explicit source term, such as in 
the single-particle density matrix, ^^-s^iiiii^ii^i^ Wigner 
functioii^^^'23^'25^^^^ and Pauli equatioii^o^ 
transport formalisms, or via a special self-energy 



term in the ubiquitous nonequilibrium Green's function 
formalism.— In this work, the problem of 
contact-induced decoherence is treated using the open 
systems formalism:^^ we start with a model interaction 
Hamiltonian that describes the injection of carriers from 
the contacts, and then deduce the resulting nonunitary 
evolution of the active region's many-body reduced sta- 
tistical operator in the Markovian approximation. The 
following two features distinguish this paper from other 
recent works r^iiSiiid^ in which Markovian rate equations 
have also been derived for tunneling nanostructures: 

(1) Derivation of the Markovian evolution is achieved 
by coarse graining of the exact short-time dynamics in 
the presence of memoryless contacts, rather than utiliz- 
ing the weak-coupling and van Hove limit,— or the 
high-bias limit ^22, Namely, electron-electron interaction 
is typically the leading inelastic scattering mechanism in 
the contacts. If the contacts' energy-relaxation time r 
due to electron-electron scattering is sufficiently short, 
then on the timescales coarsened over r, the contacts ap- 
pear memoryless and the evolution of the current-limiting 
active region can be considered Markovian.'*'^''*'* The ap- 
proximation of a memoryless environment, as applied to 
nanostructures, will be discussed in detail in Sections IIIII 
and[IVl 

(2) A model contact-active region interaction is intro- 
duced to describe the injection of carriers through the 
open boundaries and supplant the resonant level model. 
Namely, for tunneling nanostructures, like a resonant- 
tunneling diode, it is common to adopt the resonant- 
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level models^ when trying to separate the active region 
from the contacts: the active region is treated as a sys- 
tem with one or several discrete resonances. But the 
resonant-level model for the active region is inapplicable 
away from the resonances, and cannot, for instance, cap- 
ture the current increase in a resonant-tunneling diode at 
high biases (larger than the valley bias) that is due to the 
continuum states. Also, it is not a good model for simple 
structures without resonances, such as an nin diode or 
the channel of a MOSFET. So we introduce an alterna- 
tive model Hamiltonian that does not assume resonances 
a priori exist and that works both near and far from res- 
onances. It captures the open boundaries and naturally 
continuous spectrum of a nanostructure's active region, 
and describes carrier injection in a manner conceptually 
similar to the explicit source terms in the single-particle 
density matrix or Wigner function techniques. 

The paper is organized as follows: in Sec. |TT1 we 
overview the basics of the partial-trace-free formalism'*^ 
for the treatment of open systems pi Ap and present 
the main steps in the derivation of the non-Markovian 
equations with memory dressing pi B^ .^^ In Sec. IIIIl 
we discuss how the fast memory loss due to electron- 
electron scattering in the contacts can be used to justify 
a Markovian approximation to the exact evolution of the 
active region in a small semiconductor device or a ballis- 
tic nanostructure. In Sec. IIII Al we then perform coarse- 
graining of the exact non-Markovian short-time dynamics 
of an abstract open system (details of the derivation of 
the exact short-time dynamics are given in Appendix A) 
over the memory loss time of the environment in order to 
obtain a Markovian map, and we discuss the necessary 
conditions for this procedure to hold. In Sec. IIVI we 
introduce a model contact-active region interaction ap- 
plicable to a generic two-terminal nanostructure, which 
describes carrier injection from the contacts. This model 
interaction does not require that the structure a priori 
possesses resonances. In Sec. IIV Al wc formalize the re- 
quirements for the current-carrying contacts to be consid- 
ered a memory less environment. Starting from the model 
interaction and using the Markovian dynamics derived, 
we then proceed to derive the Markovian evolution and 
steady-state values for the distribution functions of the 
forward and backward propagating states in the active 
region of a nanostructure, and we also give the result 
for the steady-state current (Sec. IIVB|) . We discuss the 
relationship of the presented approach to the Landauer- 
Biittiker formalismi^iiai^i^ in Sec. HVCl In Sec. HVDl 
we work out the example of a one-dimensional double- 
barrier tunneling structure. The nonequilibrium steady 
states obtained as a result of the Markovian evolution 
at different biases produce an I-V curve that shows all 
the prominent resonant features, and we compare the 
results to those predicted by the Landauer-Biittiker for- 
malism. The manuscript is concluded in Sec. |Vl with a 
brief summary and some final remarks on the inclusion 
of scattering and lifting the Markovian approximation. 



II. THE FORMALISM 

A. Decomposition of the Liouville space 

Let us consider an open system S, coupled with the 
environment E, so that the composite SE is closed. For 
a ballistic nanostructure, S would represent the active 
region, while E would be the contacts; more generally, 
if scattering due to phonons occurs within the active re- 
gion, phonons should also be included as part of EM' 
S, E, and SE are assumed to have finite-dimensional 
Hilbert spaces, of dimensions c?s, ds, and dsdE, respec- 
tively. Consequently, their Liouville spaces - the spaces 
of operators acting on the above Hilbert spaces - are of 
dimensions dg, d\, and d'^d^, respectively. The total 
SE Hamiltonian Ti. is generally a sum of a system part 
Is ® '^5, an environment part He i8) Is, and an interac- 
tion part Tiint- The total Hamiltonian Ti. (acting on the 
SE Hilbert space) induces the total SE Liouvillian £ 
(acting on the SE Liouville space) through the commu- 
tator, which governs the evolution of the SE statistical 
operator p according to the Liouville equation 

^ = -i[n,p] = -iCp. (1) 

Ti. and C are given in the units of frequency. Dynamics 
of the open system S is described by its reduced sta- 
tistical operator ps, obtained from p by tracing out the 
environment states 

PS = Tr^p. (2) 

In general, the dynamics of ps is not unitary. A common 
approach to calculating the evolution of ps is by using 
projection operators^Si^iS that act on the SE Liouville 
space. Typically, an environmental statistical operator 
is chosen to induce a projection operator P by Pp — 
Pe ® Tre/i, where p is any vector from the SE Liouville 
space. Coupled equations of motion for Pp and (1 — P)p 
are then solved, often in the weak-coupling limit, and 
the reduced dynamics is obtained from ps — Tt^eP = 
TiEiPp). 

Most often, the projection operator utilized is in- 
duced by the initial environmental statistical operator 
p£;(0)— . The reason is that, in the most common ap- 
proximation of initially decoupled S and E, described by 
p(0) = Pb(0) (g) ps{0), the projection operator induced 
by Pe (O) will eliminate a certain memory term occurring 
in the evolution of ps- However, the result for the final 
dynamics must not depend on the projection operator 
used, as projection operators are, after all, only auxil- 
iary quantities. In this paper, we will follow the work on 
the partial-trace-free approach of Ref. i46i, that uses the 
projection operator P induced by the uniform environ- 
ment statistical operator Pe = (i^^diag(l ... 1). P has a 
unique property: it is the only projection operator that 
has an orthonormal eigenbasis in which it is represented 
by a diagonal form. Its unit eigenspace, of dimension 
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FIG. 1: (Color online) Decomposition of the total SE Li- 
ouville space into the eigenspaces of the projection operator 
P, induced by the uniform statistical operator Ji^. The unit 
eigenspace is equivalent to the system Liouville space, where 
the equivalence is given by the isomorphism (|3}. Working 
within the eigenspaces of P removes the need for performing 
the partial trace over the environmental states. 



(ig, is a mirror-image of the Liouville space of the open 
system S. Projecting onto the unit eigenspace of P is 
equivalent to taking the partial trace with respect to en- 
vironmental states^, because for any element of the SE 
Liouville space it holds 



,-1/2 



(Tr_E/i) 



a/3 



(3) 



Here, the unit-eigenspace of P is spanned by a basis \a(3) , 
while the Liouville space of S is spanned by where 
the two bases are isomorphic through the following simple 
relationship 



^\ - ^-1/2' 



(4) 



\ia,jP) is a basis in the SE Liouville space, induced by 
the bases \ij) and \aP) in the environment and system 
Liouville spaces, respectively. 

Decomposition of the SE Liouville space into the two 
eigenspaces of P (depicted in Fig. [ij is the essence of 
the PTF approach: every vector fi from the SE Liouville 
space can be written as a column /i = [^i /i2]'^, where 
/ii belongs to the unit eigenspace of P and represents 
(up to a multiplicative constant ^/d^) the system's re- 
duced component of /i, i.e., = Tr^/i = ^JdEpii- The 



other component, /i2, belongs to subspace 2 (the zero- 
eigenspace of P), where the correlations between S and 
E reside. It is important to note that the elements of 
subspace 2 (blue subspace in Fig. [T|) have zero trace over 
environmental states. 

In a similar fashion, an operator A acting on the SE 
Liouville space has a block- form with submatrices Apq^ 
p,q = 1,2, where An would be the system's reduced 
component of this operator. For instance, the block form 
of the SE Liouvillian L is given by 



L = 



•^11 ^^12 
-^21 ^22 



(5) 



where Ln is commutator-generated, and corresponds to 
an effective system Hamiltonian Hs-l-Tr£;(Hint)/rf£;- Off- 
diagonal, non-square Liouvillian submatrices, £12 and 
C21 = ^^12 1 represent the S-E interaction as seen in the 
composite Liouville space - when Tiint vanishes, so do 
£12 and £21- £22 can be perceived as governing the evo- 
lution of entangled SE states, and tends to a form fixed 
by TLs and T-Le when the interaction is turned off. 



B. Equations with memory dressing 

Using the notation introduced above, the evolution of 
the reduced statistical operator ps can be represented by 



Ps(t)^Uii{t,0)ps{Q) 



iUi2{t,0)p2iO), (6) 



where Un and U12 are the submatrices od the SE evolu- 
tion operator U, given by 



l((t,0) 



exp 



£11 £12 
£21 £22 



dt 



Un{t,0) Ui2{t,0) 

U2l{t,0) U22{t,Q) 



(7) 



In Ref. |4] 

derived as 



dUii 

dt 
dUi2 

dt 



equations of motion for Un and U12 were 

-i(£ii -£l27^)^^ll, (8a) 
-^(£ll-£l27^)^^l2-^£l2V, (8b) 



accompanied by the initial conditions Z^ii(0, 0) = 1 and 
^12(0, 0) = 0. Quantity TZ is the so-called memory dress- 
ing, as it appears to " dress" the real physical interaction 
£12 and yield an effective (generally complex) interac- 
tion term, — £i27?., that accompanies the hermitian term 
£11, responsible for unitary evolution. Memory dress- 
ing describes the cumulative effect of the S — E inter- 
action, as witnessed by a quadratic feedback term in its 
self-contained matrix Riccati^i^ equation of motion (be- 
low). The other new quantity occurring in ([8]), V(t, 0), 
can be perceived as the evolution operator for the states 
from subspace 2, and is important for the description of 
the influx of information from E to S. TZ and V obey 
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dR 

~dt 
dV 
'dt 



= —i£22'R- — iTlCuTZ + iTZLii + iC2i, (9a) 



(9b) 



accompanied by 7^(0) = and V(0, 0) = 1. 

Equations ^ and ([9]) are exact they are an alterna- 
tive form of the SE Liouville equation ((T|) . The resulting 
exact evolution of the reduced statistical operator can be 
expressed through the following differential equation of 
motion 



dpsjt) 
dt 




Ci2'R{t)] ps{t) 
dEV{t,0)p2iO). 



(10) 



which is a partial-trace-free form of = Ti E{—iLp). 

If we restrict our attention to the evolution starting 
from an initially uncorrelated state of the form 



Pi0) = pE{0)®psiO), 



(11) 



it is possible to completely reduce the problem to sub- 
space 1. Namely, it is possible to write 



P2(0) = Mpi(0)=d^'/'Mps(0), 



(12) 



where the mapping A4 is completely determined by the 
components of pe{0), the initial environment statistical 
operator (see Appendix [A|. Equation (fT2|) embodies the 
argument made by Lindblad^® that a subdynamics ex- 
ists only for an uncorrelated initial state, because, as a 
consequence of ([6]) and (fT2|) . it is possible to write 



psit) = [Uii{t,0)+Ui2{t,0)M]psiO), 



(13) 



W(t,0) 



so the evolution is completely described on the Liouville 
space of the open system. When (fT2|) is substituted into 
(fTU)l . we obtain the differential form of as 



dpsjt) 
dt 



-i [Cn - CuTZit)] ps{t) - iCi2Vit)Mps{0). 

(14) 

It is well known that a subdynamics can also be obtained 
for the case of an initially decoupled SE state by sim- 
ply choosing the initial environmental statistical opera- 
tor pe sls the one to induce the projection operator P 
(see, for instance, Ref. fss!). However, the result for the 
final dynamics must not depend on the projection op- 
erator used, as projection operators are, after all, only 
auxiliary quantities. While the physics must be the same 
regardless of the projection operator used, the opacity of 
the equations obtained certainly varies. Equation 
shows explicitly how the subdynamics looks for P; by 
generalizing the proof in Appendix |Al one can write the 
subdynamics for any other projection operator instead of 
P. The reason we are using P instead of the projection 



operator P induced by the initial environmental statisti- 
cal operator is that, as stated previously, P is the only 
projection operator that has an orthonormal eigenbasis 
in which it is represented by a diagonal form While 
any other projection operator P still projects onto its 
own (i|.-dimensional image space (see Appendix |^ , P 
and 1 — P never assume simple diagonal forms, so after 
projecting one still needs to explicitly take the partial 
trace, which leaves the equations less transparent. 



III. DECOHERENCE IN THE PRESENCE OF A 
"MEMORYLESS" ENVIRONMENT 

The non-Markovian map yV(i, 0) = Uii{t,0) + 
Ui2{t,0)M that defines the subdynamics (fT3|) can quite 
generally be written as 



W{t, 0) = T'^ exp / JC{t') dt' 



(15) 



Here, /C(t) is the generator of W{t,0). In general, 
IC(t) = —iCcS — Git), i.e., it contains an effective system 
Liouvillian CcS and a correction G due to the system- 
environment interaction, which describes decoherence. In 
case of Markovian evolution, JC = —iCcS — G = const., 
and G must have the well-known Lindblad dissipator 
formSi^ in order for the map to remain completely 
positivei^i^ 

In general, it is impossible to obtain W(t, 0) exactly. If 
one is interested in retaining the non-Markovian nature of 
psp . typically an expansion up to the second or fourth 
order in the interaction is undertakeni^^ On the other 
hand, a Markovian approximation to the exact dynamics 
can be obtained in the weak-coupling and van Hove limit, 
as first shown by Daviesi^ Although the weak coupling 
limit has been used previously by several authors^^'^^ to 
derive Markovian rate equations for tunneling structures 
in the resonant-level model, this approximation is not 
generally applicable for nanostructures.^^ 

The point we wish to make here is that the Markovian 
approximation to the long-time evolution of nanostruc- 
tures can be justified more broadly, by employing the 
approximation of a memoryless environment for the con- 
tacts. Consider first the active region of a small semi- 
conductor device; a good example is the state-of-the art 
MOSFET with 45 nm lithographic gate length (physical 
gate length is estimated to be around 20 nm), found in 
Intel's 2008 Penryn processorsi^ Semiconductor devices 
are generally required to operate at (or at least near) 
room temperature, where phonons are abundant. How- 
ever, due to the active region's minuscule dimensions, 
scattering happens infrequently, so the active region does 
feature quasiballistic transport, where scattering can be 
added as a perturbation to the ballistic solution. The 
bulk-like contacts of semiconductor devices are typically 
heavily doped (e.g., ~ 10^^ — 10^" cm""^ in silicon), 
and at room temperature all the dopants are ionized; at 
such high doping densities, electron-electron scattering 
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dominates over phonon scattering as the leading energy- 
relaxation mechanism (e.g., relaxation time for electron- 
electron scattering in bulk GaAs at 10^^ cm~^ and room 
temperature is 10 fs,^^ whereas it is about 150 fs for 
polar optical phonon scattering*'^). Basically, electron- 
electron scattering in the highly doped contacts of semi- 
conductor devices ensures that the carrier distribution 
snaps into a distribution that can be considered a dis- 
placed (also known as drifted) Fermi-Dirac distribution^ 
(see also Sec. IIV Ap within the energy-relaxation time 
r « 10"'^ — 10^ femtoseconds^i^ (the actual value de- 
pends on the doping density and temperature). This 
time is very short with respect to the typical response 
times of these devices, which is on the timescales of 
TAR ~ 1 — 10 ps ("AR" stands for the active region). 
Therefore, for small semiconductor devices, on timescales 
coarsened over the energy relaxation time t of the con- 
tacts, the contact distribution function responds virtu- 
ally instantaneously, and the contacts can be considered 
memoryless, while the relaxation of the whole structure 
happens on timescales a few orders of magnitude longer. 
(A memoryless approximation must be applied with care 
to current-carrying contacts, as we will discuss in detail 
in HvH and UFCl ) 

For low-dimensional nanostructures, fabricated on a 
high-mobility two-dimensional electron gas (2DEG) and 
operating at low temperatures, the energy relaxation in 
the contacts is also governed by the inelastic electron- 
electron scatteringjSIiSS, because the phonons are frozen 
(although there are indications that acoustic phonon 
scattering may be important down to about 4 K®^). The 
near-equilibrium energy relaxation times in these con- 
tacts are much longer than in devices, falling in the wide 
range of 10" — 10"^ pg ^es, 69^70, 7i depending on the contact 
dimensionality {1DI2jIM^ or 2D^i^), carrier density, 
and temperature. Excitations with energies higher than 
ksT, such as when bias V > ksT/q is applied across 
the nanostructure {q is the electron charge), relax more 
rapidly,^^'^^ which is of particular importance in the col- 
lector contact. In low-dimensional nanostructures, there 
are also experimental indications that coupling of the ac- 
tive region to the contacts governs its evolutioni^ As 
for the typical response times of nanostructures, recent 
experiment by Naser et alJ'^ demonstrated Markovian 
relaxation in quantum point contacts on tar ~ 50 ns 
timescales at 4 K, so the ratio t/tar is still less than 
unity, but not as small as in devices. 

Still, there is enough rationale to further explore a 
nanostructure's dynamics within the approximation of 
memoryless contacts, with the understanding that this 
approximation must generally be qualified, especially for 
nanostructures at very low temperatures. We will there- 
fore proceed with deriving the Markovian approximation 
to the exact non-Markovian equation (|14p in the pres- 
ence of an environment that loses memory on a timescale 
r, presumed much shorter that the response time of the 
open system, and we will derive the relationships that the 
coarse graining time r must satisfy for the approximation 



to be consistent. Then, in Sec. IIVI we will see what type 
of constraint that puts on our energy relaxation time in 
the contacts. 

Before proceeding with the formal development, it is 
worth stressing that the importance of a Markovian ap- 
proximation to the exact evolution is great, because with 
both nanoscale semiconductor devices used for digital ap- 
plications and with DC experiments on nanostructures, 
one is primarily interested in the steady state that the 
structure reaches upon the application of a DC bias. In 
these situations, it is sufficient to employ the Markovian 
approximation to the evolution (if warranted), as it is 
correct on long timescales and will result in the correct 
steady state. 



A. Markovian evolution by coarse graining 

To practically obtain the Markovian approximation 
due to an environment that loses memory after a time 
T, we use the coarse-graining procedure: we can parti- 
tion the time axis into intervals of length r, i„ = nr, so 
the environment interacts with the system in exactly the 



same way during each interval i„ 



PS.n+l — PS,n -p 



78 



SO 



(16) 



where K,^ 



is the averaged 



value of the map's generator over any interval [i„,t„_|_i] 
[K. is reset at each i„). If the timescales are coarsened 
over r, then the term on the left of (|16p approximates 
the first derivative at t„, so the system's evolution can 
be described by 



— = K.rPs{t). 



(17) 



The above map is completely positive and Markovian 
(coarse graining preserves complete positivit3*2§.) , but still 
has little practical value, because extracting K, explicitly 
from first principles is difficult. However, if the coarse- 
graining time T is short enough, then the short-time ex- 
pansion of K, can be used to perform the coarse-graining. 
Up to the second order in time (details of the short-time 
expansion can be found in Appendix [B)) . 



K.{t) = -iC^s - 2At + o{t^), 



(18) 



where C,s = [Hs + (Hi„t),...] = Cs + [(Hi„t),...] 
is an effective system Liouvillian, containing the 
noninteracting-system Liouvillian £5 and a correction 
due to the interaction [(...) = Ti-e[pe{0) ■ ■ ■] denotes the 
partial average with respect to the initial environmental 
state /9£;(0)]. The matrix elements of superoperator A, 
in a basis a/3 in the system's Liouville space (Liouville 
space is basically a tensor square of the Hilbert space), 
are determined from the matrix elements of the interac- 



6 



tion Hamiltonian: 

.a/3 _l(/n,2 \" r/3' /a,2 \0'i:a 
^^a'fj' — 2 \\'^int/a' " l3 + \'T-int/;3' "a' 

- ((^int)^):, < + 2(Hint)S,(7^int)J " ((Hi„t)^) J - 

where p-^ are the eigenvalues of the initial environment 
statistical operator pe{0)- A has implicitly been defined 
previousljiS in the interaction picture and with the as- 
sumption of (Tiint) = 0. Here, we work in the Schrodinger 
picture and generally need to retain (Tiint) 0, which is 
important for the inclusion of carrier-carrier interaction 
in nanostructures. A contains essential information on 
the directions of coherence loss. 

If the coarse-graining time r is short enough that it 
holds 

l|A|lr«|l/;eff|l, (20) 

then the short-time expansion of JC can be used for 
coarse-graining, and we obtain 

Kr = -tC,s - At, (21) 

leading to the Markovian equation 

^Pjp. = {-^C,s^Ar)ps(t), (22) 

which is the central equation of this paper. For the 
Markovian approximation to be consistent, the sys- 
tem's relaxation (occurring on timescales no shorter than 
1/||A||t) must be much slower than the environment's re- 
laxation (occurring over t), therefore we must have 

||A||r2«l. (23) 

Conditions (I^Dl) and (I^S)) can compactly be written as 

||A||t2 «min{l,||/:eff||T}. (24) 

B. Some general considerations regarding the use 
ofEq. ((22]) 

Before we proceed to treating a concrete nanostructure 
as an example, there are several general features regard- 
ing the use of Eq. ([22]) that can be used more broadly 
than in the treatment of nanostructures. (The reader in- 
terested exclusively in decoherence in nanostructures can 
skip the rest of this section and go directly to Sec. II Vl ) 

1. Decoherence- free evolution in the zero-eigenspace of A 

Let us assume for a moment that Ccs and A commute 
(we will see two cases of this situation in Appendix[C]). If 
so, the components of ps belonging to the null-space of A 



will not decohere - they will continue to evolve unitarily, 
as the null-space of A will be invariant under Ccs- Com- 
ponents of PS corresponding to the nonzero eigenvalues 
of A will decohere until they drop to zero. So in the case 
of commuting Ccs and A, null-space of A is decoherence- 
free. For non-commuting CcS and A, this statement can 
be generalized to 

Theorem 1. // a subspace of J\f{A), the null-space 
of operator A, is also an invariant subspace of CcS, then 
it supports decoherence- free (unitary) evolution according 
to the map [2^1 . 

Proof Let N'{A) be a subspace of N{A). If TV' (A) 
is an invariant subspace of £eff) then it is an invariant 
subspace of the full generator of the Markovian semi- 
group (|22p . and consequently an invariant subspace of 
the semigroup. A statistical operator p*^, initially pre- 
pared in A/"' (A), would remain in Af'{A) at all times, and 

evolve unitarily according to ^ = — i£^gp°(t), where 
£°ff is the reduced form of £eff onto N'{A). ■ 

This theorem is equivalent to the statements made 
in the original works on decoherence- free subspaces^", 
where a decoherence-free statistical operator was defined 
through annulment by the Lindblad dissipator. Note, 
however, that here wc identify the decoherence-free sub- 
spaces in the system Liouville space, rather than in its 
Hilbert space. This allows for the possibility that some 
entangled system states (Trp| ^ Tr/95) could be resilient 
against decoherence, which is a potentially useful feature 
that cannot be captured in the Hilbert space alone. 

Theorem 1 gives us a straightforward, general recipe 
for the classification of the decoherence-free subspaces in 
the case of Markovian dynamics (P^ . What one needs 
to do is to construct the operator A according to Eq. 
(fT9)l . from the microscopic interaction Hamiltonian and 
the environmental preparation, solve its eigenproblem (in 
general numerically), and investigate whether any of its 
null-spaces is invariant under >Coff- This is a simple, ef- 
ficient way to approximately determine where the infor- 
mation should be stored, and should work well as long as 
the system is small enough to allow for a full solution to 
the eigenproblem of A. 

Moreover, the structure of the eigenspaces of A en- 
ables us to determine the directions of decoherence. For 
instance, regardless of the value of r, we can still tell 
which states do and which do not decohere, and calculate 
the relative values of the decoherence rates for two given 
states. For fast switching in nanoscale semiconductor de- 
vices, for example, we need rapid coherence loss between 
the active region and leads, and we may therefore opt to 
prepare the system in the subspace of A corresponding 
to one of its largest eigenvalues. 

2. Identification of the steady state 

An important special case of a decoherence-free sub- 
space is that of a vector belonging to the intersection of 
M{C^s) and M{A). 
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Theorem 2. A statistical operator belonging to 
Af{Ccs) n the intersection of the null-spaces 

Af{Cos) o,nd J\f{A), is a steady state for the evolution 
according to the map i22\). 

Proof M{Cos) n Af{A) is the null space of the Marko- 
vian semigroup generator. Consequently, any statistical 
operator prepared in A/'(£eff ) H A/'(A) remains unchanged 
at all times, satisfying the definition of a steady state. ■ 

By looking into the common null-subspace of both £off 
and A, one can narrow down the set of potential steady 
states, which is important in many-body transport calcu- 
lations. In the case of a many-particle open system, a full 
solution to the eigenproblem of A may not be tractable; 
however, identification of the common null-space of CcS 
and A may be. 



3. A comment on the validity of Eq. I22\l 

In general, whenever an efficient resetting mechanism 
can be defined for the environment, so that Eq. (|24p is 
satisfied, ((22l) should be applicable. However, it also ap- 
pears that the simple equation (|22p may be used more 
broadly than specified by ([24|) . Namely, on one of the 
few exactly solvable systems, the spin boson model with 
pure dephasing, which experiences Markovian evolution 
in the long time limit regardless of the coupling strength, 
it can be shown (see Appendix I C ip that one can define a 
mathematical coarse graining time r that is shorter than 
any other timescale in the coupled system and environ- 
ment, so that coarse-grained evolution over r (f22| and 
the exact Markovian evolution coincide. So, it appears 
that not only does coarse graining result in Markovian 
maps, but the converse might also be true: it is possi- 
ble that a given Markovian evolution can be obtained by 
coarse graining of the short-time dynamics if a suitable 
(ultrashort) mathematical coarse graining time is cho- 
sen. This statement would, of course, be very difficult to 
prove in general terms, but is interesting because it would 
mean that all one needs to deduce the steady state for 
the evolution of an open system is the information on 
its short-time dynamics (|18p . which can in principle be 
done relatively straightforwardly and from first princi- 
ples (the microscopic interaction and the preparation of 
the environment). Indeed, on an additional example of 
the Jaynes-Cummings Hamiltonian in the rotating wave 
approximation, which has been worked out in Appendix 
IC 2[ it has been shown that by using map and the 
resulting criterion for the steady state (Theorem 2), re- 
laxation towards the proper equilibrium state has been 
obtained. So it appears that the applicability of Eq. (1^^ 
may extend beyond the formal range of its validity (|24p . 



IV. A TWO-TERMINAL BALLISTIC 
NANOSTRUCTURE 

In this Section, we consider a generic two-terminal 
nanostructure under bias, and introduce a model 
interaction between the ballistic active region and 
the contacts. This model should hold regardless 
of whether the structure has resonances or not, as 
it is constructed to mimic the source term in the 
single-particle density matri x^^i^^i^''''^^i^^i^° and Wigner 
function^^^^S^^^^^^ formalisms, and pre- 
serve the continuity of current, state-by-state. In Sec. 
IIVDI the results are illustrated on a one-dimensional two- 
barrier tunneling structure. 

The left contact is the injector (source), biased neg- 
atively, while the right contact is the collector (drain). 
The contact-active region boundaries are at (left) and 
xr (right), with W = xr — xl being the active region 
width. We will assume that the active region includes a 
large enough portion of the contacts (i.e., exceeding sev- 
eral Debye lengths) so that there is no doubt about the 
flat-band condition in the contacts. Also, W should be 
large enough to reasonably ensure a quasicontinuum of 
wavevectors (Afc — 2tt/W) following the periodic bound- 
ary conditions. While sweeping the negative bias on the 
injector contact, we will assume that it is done slowly (so 
that between two bias points the system is allowed to re- 
lax) and in small increments (so that the potential profile 
inside the active region does not change much between 
two bias points, and can be regarded constant during 
each transient). 

For every energy £k above the bottom of the left 
contact, the active region's single particle Hamiltonian 
has two eigenfunctions, a forward (5'fe) and a backward 
(^'_fc) propagating state, that can be found by (in gen- 
eral numerically) solving the single-particle Schrodinger 
equation for a given potential profile in the active re- 
gion. To keep the discussion as general as possible, 
we will not specify the details of how the active re- 
gion actually looks (Fig. ^ - e.g., it can be a het- 
erostructure, a pn homojunction, or a MOSFET chan- 
nel - but we will require that the contact-active region 
open boundaries (at xl and xu) are far enough from any 
junctions in the active region, so that the behavior of 
near the junctions is already plane- wave like, i.e., 
that their general form near the injector is '^k{x^) = 
^ikx+ _|_ ^_^^^^-ikxl ^ ^i_k{xj^) = t^k.Le^'''^^ , while 
near the collector "iikix^^) — tki^RC^'^^R, 5'_fc(a;^) = 
g-ifc _|_ ff^, j^Q^^ Here, where i's and r's are the 
transmission and reflection amplitudes, while k and k' 
are the wavevectors that correspond to the same energy 
£k, measured with respect to the conduction band bot- 
toms in the left and right contacts, respectively (fc^ = 
k'^ — 2mqV/h^ — 2m£k/h'^, where —V has been applied 
to the left contact, and q is the electron charge). 

Associated with 5'^ (vlZ-fc) in the active region are the 
creation and destruction operators dj. and dk (c^L^ and 
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FIG. 2: (Color online) Schematic of the coupling between 
the active region of a generic two-terminal nanostructure and 
the contacts. In case of ballistic injection through the open 
boundaries, a forward-propagating state "i/k is coupled with 
the states exp{±ikx) in the left contact (injected and reflected 
wave) and the state exp{ik'x) in the right contact (transmit- 
ted wave) via a hopping model interaction 



d-k), SO the active region many-body Hamiltonian is 

Hs^Y. M4dk + dlf^d^k). (25) 

k>0 

Spin is disregarded, and tuk = Sk/h. In case of ballis- 
tic injection through the open boundaries, each state "^k 
is naturally coupled with the states exp{±ikx) in the 
left contact (injected and reflected wave) and the state 
exp{ik'x) in the right contact (transmitted wave). For 
^_fei the coupling is between exp{zLik' x) in the right and 
exp{—ikx) in the left contact. To model this coupling via 
a hopping-type interaction, we can write quite generally 
(see Fig. E]) 



k>0 



Akdlck,L 



■Aldi^^^dk + Alclj^dk 



(26) 



(c±fc,L) and c\, ^ {ck'.R) create (destroy) an elec- 



tron with a wavevector ±k in the left and k' in the right 



contact, respectively. The hopping coefficients 
and A|, are the rates of injection, reflection, and trans- 
mission, respectively. Therefore, they are proportional 
to the injected, reflected, and transmitted current for the 
state i.e.. 



A[, A* ^ 



(27) 



where TLk and 7^ are the reflection and transmission co- 
efficient at a given energy. The actual magnitude of Afc 
can be determined by requiring that A| = T^A^, the 
hopping rate from the active region into the right con- 
tact, be the same as the current (per unit charge) carried 



through the active region by ^fc. This just means there 
is no more reflection once the wave exits the active re- 
gion and gets into the outgoing contact, and is usually 
referred to as the assumption of reflectionless leads. 
The current carried by '^k is given by the well-known 
quantum-mechanical relationship 



Ik 



qh 



qh\tk',R\^k' 



^Imin^^k) 

qhk Tk 



\^k\ 



(28) 



where we have used the form of ^'fc near the right con- 



tact tk' .Re"^' , and ||^'fe||^ = dx\^k{x)\'^ is the norm 
squared of "^k over the active region of width W . Since 
we require that Ik/q = A^ = T^Afc, we find 



-w 



Ak 



hk 



mll^fell 



(29) 



Finally, 

(30) 



k>0 



while the Hamiltonian for the backward propagating 
states can be written in an analogous fashion, as 



with A_ 



A-k (dlc^k',R + TZ-kc\, j^dk + T^kc\, j^di 



k>0 



' h.c. 
(31) 



nk' 



and % =T-k. T^k 



7^_ 



m||*-fclP ' 

When we put it all together, we have for the interac- 
tion Hamiltonian of the active region with the left/right 
contact: 



h.c. ntt = 5]Afc{(c^ + 7^fccL,,^)4 



fc>0 

4 {ck,L + TlkC-kx 



)} 



(32a) 



k.Ld-k+'T'kdl^c_k,L 



} 



n_/R 
'T-int 



J2^~k{{clk',R + T^kcl^i,) d_ 

fc>0 

{C-k'M. + Ti-kCk'M)^ 

Ak ^Tkc\, j^dk + Tkd\ck',R^ . 



(32b) 



A. Current-Carrying Contacts and the 
Approximation of a Memoryless Environment 

Now that we have the interaction Hamiltonians in 
place, we should evaluate the matrix elements of the su- 
peroperator A, which leads us to the questions how the 
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approximation of a memoryless environment is actually 
applied to contacts carrying current, and how the expec- 
tation values of the interaction Hamiltonian are to be 
calculated. 

In general, as the current flows through the structure, 
we must allow for different distributions of the forward 
and backward propagating waves in the left and right 
contacts to ensure current continuity. A simple and often 
employed approximation for the steady-state distribution 
in the contacts carrying current / is a single-parameter 
drifted (or displaced) Fermi-Dirac distribution^Sii^ 

f^k = {n±k,L) = ■ (33) 

Here, kp is the Fermi wave vector and kd is the drift 
wavevector, determined from the total current / as 
kd = ml /riioqh, where nuj is the ID carrier density in 
each contact (contacts are assumed identical). A drifted 
Fermi-Dirac distribution, with the temperature equal to 
that of the lattice, is often employed when we are in- 
terested in just the first two moments of the distribu- 
tion function (i.e., maintaining charge neutrality and en- 
suring current continuity). Additionally, if needed, in- 
formation on the electron heating can be incorporated 
in this distribution by allowing for a discrepancy be- 
tween the electronic and lattice temperatures (we will 
neglect electron heating here) . Detailed ensemble Monte 
Carlo - molecular dynamics simulations of carrier trans- 
port in highly doped (> 10^'' cm'^) bulk semiconductors, 
in which electron-electron scattering is the most efficient 
energy relaxation mechanism, have shown to produce dis- 
tributions very close to the drifted Fermi-Dirac distribu- 
tion (pg|) |6ii63,82 .y^rj^ich is generally accepted as a decent 
approximation for these systems. Here, we will also adopt 
(|33p for the distribution of carriers in the current-carrying 
contacts, and it is reasonable if the (one-dimensional) 
contacts are longer than V Dt, where D is the diffusion 
constant (otherwise, the distribution function in them 
may not be thermalized^i^i) . 

Now, the question arises what happens if we try to 
sweep the voltage. We have mentioned before that the 
voltage is to be swept slowly (enough time between two 
bias points for the system to relax) and in small incre- 
ments (so that we can consider the barrier as having a 
constant profile during each transient). The latter is cru- 
cial for the implementation of the approximation of a 
memoryless environment. Suppose that, at a bias V, a 
steady-state current / is flowing through the structure. 
If we increase the bias to AV a.t t — 0, where AV is 
very small, within the first t — t, the current is vir- 
tually unchanged - it takes the current a much longer 
time TAR ^ r to change significantly ( AR stands for " ac- 
tive region" ; once we have had a chance to complete the 
calculation, we will see that tar will be equal to l/Ar, 
where A is a relevant eigenvalue of A). Therefore, after r, 
the contact carriers have redistributed themselves to the 
old distribution function that they had at V. Basically, 



the contact carriers as redistributing themselves to Eq. 
P3p determined by the (virtually) instantaneous current 
level at each r; the current, however, changes very little 
during each t. By the time the current has saturated 
(~ Tar) J the contact carriers have had a chance to get 
redistributed many times; however, if the total voltage 
increase AV is very small, the total current increase dur- 
ing the full transient will also be small, so we can say that 
during the whole transient the distribution functions of 
the forward and backward propagating states have been 
resetting to nearly the same distribution, approximately 
the average of f^'^ over the interval [V, T^-l- AT^] . Clearly, 
as the voltage sweep increment AV 0, we can say that 

^ P _____ 

during a transient the contacts redistribute to /^j, ([33]) 
at V. 

Evaluation of kd that enters the contact distribu- 
tion functions at a given voltage must be done self- 
consistently: starting with a guess for kd at a given volt- 
age, steady-state distributions and current are evaluated 
(as detailed in the next section) . The obtained current is 
then used to recalculate kd, and the process is repeated 
until a satisfactory level of convergence is achieved. (Of 
course, the initial guess for kd at any voltage can be 
kd = 0, but for faster convergence it is better to start 
with the kd found for the preceding voltage.) 



B. Markovian Relaxation for a Two- Terminal 
Nanostructure. Steady-State Distributions and 
Current 



Since the interaction Hamiltonians (j32|l arc linear in 
the contact creation and destruction operators, and we 
can approximate that each contact snaps back to a 
"drifted" grand-canonical statistical operator, we have 

L/ R 

(Hjjjt ) = 0. This means that Cs = ^oS, and also leaves 
us with only the first three terms in Eq. (fT9|) for A to 
calculate. One can show that A = A^ + A^, where 

- E^l«(^.^n{^)L"(Wfn('')f;- (34) 



10 



The first and the second term in Equation ((34)) give a where 



general contribution of the form A"^, since 



A:>0 

+ [{l-{nk,L))+nUl-{nk^L))]dldk} 
+ ^-k{TK^{n-k,L)d-kdU 
+ T,^{l-{n_k.L))dl„d_k} 
= J2Al{{ft+Tllf^l,)d,dl 

k>0 

+ [ii^f^) + nl{i-f\)]dld,} 
+ r,2(i-/^,)dV_fe} 



(35) 



preserves the filling of states. We have used {n±k^i) = 
f^f., where is the drifted Fermi-Dirac distribution 
function in the left contact (l33l). 



In contrast, the third term in (|34p 



= j2Ai{[ii- ft) + niii f\)] {d,r^, (dif; 

+ {ft+nf'^k)idl)2'idkfp} (36) 

+ A^_,{r,2(i-/f,)(d_,.)"'(rfU)f 



f^kidlX'id-kf,'} 



gives a contribution of the form A^^. 

Each term in A is a sum of independent contributions 
over individual modes [A = J2k ^k] that attack only 
single-particle states with a given k. The same holds for 
Cs- Consequently, in reality we have a multitude of two- 
level problems (see Appendices IC ll and IC 2p . one for each 
state ^k, where the two levels are a particle being in '^k 
("-I-") and a particle being absent from '^k ("-")• Each 
such 2-level problem is cast on its own 4-dimensional Li- 

ouville space, with pk = {Pk~^ ^ Pk~ ' Pk~^ ^ Pk~) being 
the reduced statistical operator that describes the occu- 
pation of According to (P^ . 



dpk 
dt 



-i^s,k - AkT]pk, 



(37) 



S,k 



Ak = 



■ 








0" 





2ujk 














-2ujk 





. 








0. 



Ak -Bk 

Cfc 

Cfc 

-Ak Bk 



(38a) 



(38b) 



and Ak = AUil - ft) + 7^2(l - /f,) + T,'{1 ft:)}. 
Bk - Al{ft+nlfb^+T^f§}, and Ck = {Ak + Bk)l2 = 



A?-(l 



TZL +Tt) /2. The rows/columns are ordered as 



1 = |+)(+|,2 = |+)(-|,3 = |-)(+|,4 = l-X- , 
The diagonal elements in A^ originate from the terms 
of the form A^^, calculated using Eq. pS)) . while the off- 
diagonal ones originate from A^^ (155)) . Strictly speak- 
ing, the time evolution above is valid if (|24)) is satisfied, 
which in this case implies A^^r^ <C min{l, w^r}. Af- 
ter approximating « W, we obtain the condition 
[vt/W)'^ ^ min{l, Wfcr}, where v = hk/m. For typ- 
ical values of = lOOnm, v < vp — 10^ rn/s, and 
m = O.Q67mo appropriate for GaAs, both equations will 
be satisfied for r ^ 1 ps. 

Clearly, off-diagonal elements pt~ and p^T"'" decay as 
exp (=Fi2tjfc — TCk)t and are zero in the steady state. The 
two equations for p^"*" = fk{t) and p'^~ = 1 — fk{t) are 
actually one and the same, and either one yields 



dfkjt) 
dt 



-TAkfk{t)+TBk[l- fk{t)] 

-T{Ak + Bk)fk{t) + TBk, 



(39) 



where fk is the distribution function for the active region. 
In the steady state, we have — ^^^'-^^ (for — fc, by 
analogy), so finally 



roc 

Jk 

rOQ 

J-k 



ft 



nfik+Ttft 



1 



Tt 



(40a) 
(40b) 



Note that there is no dependence of the steady-state 
distributions on A^^, the hopping interaction strength, 
or the coarse-graining time r. obviously differ from 
the contact distributions (sec discussion in the next sec- 
tion). The discontinuity of the distribution functions 
across each open boundary is a price to pay to conserve 
the flux across it, the same as in the heuristic treatment 
of carrier injection in the density matrix, Wigner func- 
tion, and Pauli equation formalisms (see the discussion 
on p. 4907 of Ref. [Ml). 

The steady-state current (per spin orientation) can be 
calculated as 



k>0 



f^h + f^kl- 



(41) 
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where h = qhkTk / m\\^ and I-k = qhk'Tk/m\\^- 



Ij-k are each constant across the active region and 
given by . The total current carried by the forward prop- 
agating states (per spin orientation) is 

roc _ if^sr f kTk 



fc>0 



m 2tt 



(42a) 



where we have used kdk ^ mdEk/fi^ and Afc « W/2tt. 
Similarly, the current component (per spin) carried by 
the backward propagating states is 



1° 



qW 



d£k f- 



(42b) 



so the total current (per spin orientation) can be found 
as 



outgoing states {f^^ ^'^'^ f§) contacts in order to 

calculate the distributions of the forward and backward 
propagating states {f^ and /^,) in the active region. 
The reason is that the information about the outgoing dis- 
tributions, supplied by the transfer matrix, is destroyed in 
the contacts, where the inelastic scattering very rapidly 
redistributes carriers. 

Our model for the inelastic current-carrying contacts 
can actually be considered as complementary to the well- 
known model of voltage probes i^^'^^i^^'^'^ On average, a 
voltage probe carries no current. Due to inelastic scatter- 
ing, the distribution function in a voltage probe is reset 
to the equilibrium one on timescales much shorter than 
the response time limited by the active region (tar)- In 
contrast, there is no voltage drop over a current-carrying 
contact (conduction band bottom is flat), while the aver- 
age current carried by it is generally nonzero. Due to in- 
elastic electron-electron scattering, the distribution func- 
tion in the current-carrying contact is reset to a displaced 
Fermi-Dirac distribution on timescales much shorter than 

TAR- 




foe 



(43) D. Example: A Double-Barrier Tunneling 

Structure 



This expression is parameter-free, because ||5'±fe|| in the 
denominator scale with W. 



Relationship to the Landauer-Biittiker 
formalism 



A natural question emerging at this point is how the 
current relates to that predicted by the Landauer- 
Biittiker (LB) formalism^i^'^i^ (comprehensive re- 
views of the LB formalism can be found, for instance, 
in Refs. [s^ and IsgI. as well as in many textbook o^^i^^ ). 
The one-channel variant of the current formula is referred 
to as the Landauer formula. 



TOO Jl_ 

Lan - ^ 



d£k [/^(ffc)-/"(ffc)jTfc, 



(44) 



where f^{£k) and f^{£k) are the equilibrium distribu- 
tions in the left and right reservoirs j2§^ Generalization to 
multiple channels is due to Biittiker ;^°i^^'^^i^^ 

Both the LB approach and the approach presented here 
focus on maintaining the carrier flux through the open 
boundaries between the active region and the contacts. 
There is one major difference, however. In the LB ap- 
proach, what is known are the distributions of the states 
entering the active region (in our notation, and f^^')' 
nothing is specified about the distributions of the states 
going out of the structure {f^i. and /^), as they can be 
calculated by using the transfer matrix (a nice exposi- 
tion of this issue can be found in Refs. [oO and IsgI). In 
contrast, in the approach presented here, we need the in- 
formation on both the incoming (/^ and f^^') ^^'^ 



We illustrate the results of Sec. IIVBI on a one- 
dimensional, double-barrier tunneling structure, formed 
on a quantum wire in which only one subband is popu- 
lated. The Fermi level is at 5 meV with respect to the 
subband bottom. The well width is 15 nm, the barrier 
thickness is 25 nm, and the barrier height is 15 meV. 
These result in one bound state at about 6.84 meV when 
no bias is applied. The goal is to calculate the nonequi- 
librium steady-state distribution functions specified by 
Eq. (I40p under any given bias V, and use this informa- 
tion to construct the I-V curve. For simplicity, in this 
calculation the voltage is assumed to drop linearly across 
the well and barriers, but in general, Eqs. (|40|l need to be 
coupled with a Poisson and a Schrodinger solver to obtain 
a realistic potential profile and charge distribution. 

Figure [3] shows the I-V curve of the double-barrier 
tunneling structure, as calculated according to the ex- 
pression (|43p and the Landauer formula In the 
voltage range depicted, the current flowing through the 
structure is so low that the equilibrium distribution func- 
tions in the contacts (/^(£fc) and f^{£k)) and the drifted 
Fermi-Dirac distributions p3|) . with kd determined self- 
consistently, are extremely close to one another, and give 
almost identical (|40|) and the values for the current 
The difference between the curves obtained by us- 
ing the equilibrium contact distributions and the drifted 
Fermi-Dirac is barely visible within the voltage range pre- 
sented (the maximal difference between the currents ob- 
tained these two ways is ~ 1Q~^^A). 

Both (f43|) and (f44|) describe ballistic transport, so no 
crossing of the curves typical for the inclusion of inelastic 
scattering should be expected (inelastic scattering causes 
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FIG. 3: (Color online) I-V curve for the double-barrier tun- 
neling structure, according to the expression (|43|l (solid curve) 
and the Landauer formula (|44l) (dashed curve) at 1 K. 
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FIG. 4: Steady-state distribution functions (|40p of the for- 
ward (/fe°) and backward (/^) propagating states, at the 
peak voltage from Fig. [3] (5 meV) and 1 K. Energy is mea- 
sured with respect to the bottom of the injector contact (Fig. 
[2]). Significant deviations from the Fermi-Dirac equilibrium 
distributions in the contacts coincide with the peak in trans- 
mission. 



the peak to lower and the valley to rise, so the curves 
cross^^). Both curves in Fig. [3] properly display the res- 
onant features, but the Landauer formula (|44p predicts 
a higher peak current than (|43)l . The reason is that /^,, 
used in (|43p . coincide with the contact (nearly equilib- 
rium) distribution functions only if the transmission is 
not high. Near a transmission peak, significant devia- 
tions of (|40|) from the contact distribution functions 
occur, as shown in Fig. [4] for the peak voltage from Fig. 
m and lead to the lowering of the current observed in Fig. 

m 



V. SUMMARY AND CONCLUDING REMARKS 

In this paper, a simple theoretical description of the 
contact-induced decoherence in two-terminal nanostruc- 
tures was provided within the framework of the open 
systems theory. The model active region - contact in- 
teraction was introduced to ensure proper carrier injec- 
tion from the contacts. The steady-state statistical op- 
erator of the active region was calculated by relying on 
the Markovian map derived through coarse graining of 
the exact short time dynamics over the energy relaxation 
time of the bulk-like contacts. The ballistic-limit, steady 
state distribution functions of the forward and backward 
propagating states for a generic two-terminal nanostruc- 
ture have been derived. The approach was illustrated 
on the example of a double-barrier tunneling structure, 
where an I-V curve that shows all the prominent reso- 
nant features was obtained. The relationship between 
the present approach and the Landauer-Biittiker formal- 
ism was addressed. 

The inclusion of scattering within the active region 
would alter the form of CeS, while scattering between 
the active region and the contacts (e.g. phonon-assisted 
tunneling) would essentially alter A. Equations ([40]) 
are the ballistic limit of the active region's nonequilib- 
rium steady-state distributions, and are a better start- 
ing point for transport calculations with scattering than 
the equilibrium distributions: for instance, the single- 
particle density matrix ^'■^•'(fci, ^2) = Tr sidl^dkiPs) 
in the ballistic limit is obviously diagonal, so to in- 
clude scattering within the active region, one simple 
way would be to follow the single-particle density ma- 
trix formalism, i^iilii^iiaiSa with the diagonal p^^\k, k) 
specified by (|40|) as the ballistic limit. [Clearly, p^^'i would 
no longer be diagonal in k once scattering is included.] 
Scattering due to phonons within the active region is 
generally amenable to a weak-coupling approximation, 
so it can be treated as a perturbation within the Born 
approximation. To treat phonon-assisted injection from 
the contacts, the contact many-body Hilbert space can 
be augmented to formally include a tensor product of the 
contact and the phonon Hilbert spaces , ^"'^"'^ but again a 
simpler perturbative treatment may be enough. As for 
the treatment of electron-electron scattering, A is in the 
form that allows for its inclusion between the active re- 
gion and the contacts, but this is likely to be a difficult 
technical issue. 

Finally, an important feature of the present approach is 
that it can be, at least in principle, extended to arbitrar- 
ily short timescales by forgoing the coarse-graining pro- 
cedure, so non-Markovian effects can be observed. How- 
ever, since the coarse-graining procedure phenomenolog- 
ically accounts for the efficient electron-electron interac- 
tion in the contacts, without it we would be required to 
explicitly include this interaction in the contact Hamil- 
tonian, which will require certain modifications to the 
present approach. 
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APPENDIX A: UNCORRELATED INITIAL 
STATE AND THE EXISTENCE OF A 
SUBDYNAMICS 

In this Appendix, for an uncorrelated initial state of 
the form p(0) = pe ® Ps{0), we will exphcitly show 
that /92(0), the component of p(0) belonging to the zero 
eigenspace of P, can be written in terms of ps{0) via 
equation p^ . repeated here 



P2(0) -Xpi(O) = d^5Xps(0). 



(Al) 



Together with ([6]) , this equation proves that a subdynam- 
ics exists. We will explicitly derive the mapping A4 that 
is uniquely fixed by p^;. 



1. Eigenbasis of P 



Let us first remind ourselves of the structure of the 
eigenspaces of P. Its unit eigenspace is d|-dimensional, 



spanned by vectors of the form 



kI3) = 



dm 

E 



(A2) 



This form holds regardless of the environmental basis 
chosen, which is in agreement with the fact that the uni- 
form environmental statistical operator (the one that 
induces P) is a scalar matrix, i.e., diagonal in any en- 
vironmental basis. In the zero eigenspace of P, for any 
choice of the environmental basis, we can identify two 
subspaces: 

1) A subspace spanned by vectors of the form \ia,jf3), 
with i j. This subspace is dgdEidE — l)-dimensional. 

2) A subspace spanned by linear combinations of \ia, z/3), 
which are orthogonal to all | a/3) . These are given by 



\bi,al3) 



1 



dE 



for every pair a, (3 and for i = 1, 
space is d'g{dE — l)-dimensional. 



dE + l~i j 

(A3) 

. . , c?£; — 1. This sub- 
Note how the coeffi- 
cients in the linear combinations do not depend on 

2. Range (Image Space) of P and Null Space of P^ 

Now let us get back to the initially uncorrelated state 
of the form p(0) = pe® Ps{0), and choose the eigenbasis 






system Liouville 




space 
























basis 










P=0 



SE Liouville space 
{decomposition via } 



SE Liouville space 
(decomposition via ) 



FIG. 5: (Color online) Two decompositions of the total SE 
Liouville space: the first one (right) is into the eigenspaces of 
P, induced by the uniform statistical operator pg- (System 
Liouville space and P's unit eigenspace are equivalent.) The 
other decomposition is into TZ{P), the range of the projector 
P induced by the initial environmental statistical operator 
pE, and JV{P^), the null space of P^ and orthocomplement to 
TZ{P). These two decompositions enable us to construct the 
mapping Ai used in the reduced dynamics (|14|l . 



of pe as the environmental basis — 1, . . . rfg. pE, 

the initial environmental statistical operator from p(0) = 
Pe(^Ps{0) induces its own projection operator P, so that 
for any vector p from the SE Liouville space we can write 



Pp ^ PE® (Ttep) 



i,a',/3' \ k I 



id) (A4) 



The Latin indices i^k — l,...dE count environmental 
states, while the Greek ones a,/3 = l,...ds count the 
system states. The essence of the following proof is to 
write any Pp in terms of the eigenvectors of P, and then, 
since Pp{0) = p(0), draw important conclusions about its 
components p2{0) and ps(0). 

P is not Hermitian or diagonalizable. We can, how- 
ever, still speak of its range (space of images) TZ{P), to 
which p(0) belongs becaue Pp{Q) — p(0). The orthocom- 
plement to 7^(P) is A/'(Pt), the null space of the adjoint 
operator . It is easily noted that all vectors of the form 
\ia,j(3), with i j, are in the null spaces of P, P and P^. 
Therefore, A/'(P^) is at least d^dE{dE — l)-dimensional. 
Where is the rest of A/'(P^), i.e., what is a general form 
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of a vector 

= Y.^l I*"' ' (^'^) (Cg,a/3|PM> = 0? (A5) 



indeed constitutes an orthonormal basis in Ti{P), and 
that 



(c,,a;3|PM> = £(C;)Vij(*a,z/?ba',j/3')(Tri;M)"' 



(A6) 



aP\Pfij =0 iff TiB^ = (A9a) 
ap\cg,c,'p') =0, ya,P,q,a'p'. (A9b) 



Therefore, 

{c,^.p\Pl^) = ^ Y.^^irp% = (A7) 

1=1 

Columns (f^, . . •Cg^)"'" satisfying (|A7p constitute a — 
1-dimensional space, so we conclude that A/'(P^) is of 
dimension d|d_E(rf_E — 1) + c^K^b — 1) = dgi^E ~ !)■ 
Therefore, the rank of P [dimension of Tl{P)] is so 
it is isomorphic to the unit eigenspace of P and to the 
system Liouville space. One can show that the choice 



~\ 1 - 

a/3) = -y^^Xl^B \ia,il3) 



(A8) 



-E 1=1 



Why was this analysis necessary? Because an uncorre- 
lated initial state satisfies Pp{0) = p{0), which means the 
initial statistical operator belongs completely to 7?.(P). 

Therefore, it can be written in terms of the basis a/? 
as 



(AlO) 

In Fig. [5l mutual relationships among the eigenspaces of 
P and the null and image subspaces of P are depicted. 
We obtain 



{a(3\p{0)) = (a(3\a(5) U(3\p{0)) + ^ (a/?|c,,„^) (c,,«^|p(0)) = (a(3\a(3) (a/3|p(0) 



(All) 



=0 



The important point to note is that (a/?|p(0)) and (a/3|p(0)) are equivalent up to the multiplicative constant 

aj3\aj3) 



(A12) 



We can now obtain the projection of p2(0) onto the zero-eigenspace of P as 

(&,,a;3|p(0)> - (&,,a;3M) (^|p(0)) + J] (&,a/3 |c,,a;3) (c,,a/3 1 p(0)) = (^^^|^^) (^|p(0)) 



=0 



Since ( b 



= , we arrive at 



3=1 / 



Equations above fix the mapping p2(0) = Mpi{0) from 
(jAip . and explicitly embody Lindblad's argument on the 
existence of a subdynamics^. 



APPENDIX B: SHORT-TIME DECOHERENCE 
IN NON-MARKOVIAN SYSTEMS 

In this Appendix, we formally show how to obtain the 
short-time limit to the exact completely positive non- 
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Markovian dynamical map governing the evolution of p5, 
in the form 



iat')dt' 



PsiO) 



(Bl) 



dt'[-ic,s{t')~g{t')]}ps{0), 



ps{t) - T^exp 
= T'^cxp 



where Ccg{t) is a still undetermined effective Liouvillian, 
and Q{t) is the dissipator term. It is well known that 
the form above holds for the dynamical semigroup in 
the Markov approximation, where the time-independent 
semigroup generator —iCeS — Q = const, is of the well- 
known Lindblad form^i^ that ensures the map's com- 
plete positivity. 

Here, we will perform the short-time Taylor expansion 
of the exact equation (fTl]) up to the second order in time 



o(i'), (B2) 



and we will identify the terms in the first and second 
derivatives from the desired equation (|B1[) 



dps 
dt 

d'ps 
dt^ 



-iCcs{0)-gmps{0) 
'dCcs\ f dQ 



(B3a) 



dt 



dt 



Ps(0) 



(B3b) 



with those obtained from the exact evolution described 
byEq. (fTi]). 

Indeed, by using the initial conditions i?(0) = and 
V(0,0) = 1 given in Eq. from Eq. we directly 

obtain 



(dfs\ 



V dt 



= -i{Cii+Ci2M)ps{0) 

[7^5 + (Hint), P5(0)]. 



(B4) 



Here, we have used the facts that Cn is generated by 
the Hamiltonian Hs + Tiint, where Hint = Tr£;(Hint)/rfE, 
while C12M is generated by the Hamiltonian (Hint) — 
Hint, where (Hint) — Tr^ (pfiHint) is the averaged inter- 
action Hamiltonian (see Appendix IB 1[) . Consequently, 



£eff (0) = £11 + = [Hs + (Hint),...],(B5a) 

^(0) = 0. (B5b) 

Taking the first derivative of Eq. ((T4|) and employing 



i?(0) = 0, 
^£22 [Eq. 



\ dt )q 



«£2i, V(0,0) 
we directly obtain 



1, and 



rfV(t,0) 
dt 



PS 



dt^ 



{C12C21 + C12C22M) psiO) 



- Cn{Cn+Ci2M)ps{0). (B6) 
After subtracting [-iCcsiO)]'^ psiO) from , what 

we obtain is action of the operator —i ( '^^l" ) q — 
on ps{Q)- Therefore, we will introduce operator A 



dt Jo \ dt Jo 

as 



2A 



dCcff 
dt 



d^\ 
dt), 



(B7) 



where 



2A 



{C12C21 + C12C22M) ~ C12M (£11 + £i2A^) . 



Operator A contains essential information on the di- 
rections of coherence loss in both non-Markovian and 
Markovian systems. After a straightforward calculation, 
documented in Appendix IB 11 one can obtain the matrix 
elements of A in the tensor-product basis of the system 
Liouville space 



K^P' - ^E^^|(^?nt)::<-2E(^.nt)l"(Hint)^^' + (HL)S;c| 
k { k' ) 

- ^ {((Wint)'):, < - 2(Hint)S.(Hint)J + ((Hint)') J C} , (B8) 

where, for simplicity, the environmental basis is assumed to be the eigenbasis of the environment initial statistical 
operator pE- In a more compact form, the action of A on any vector ps from the system Liouville space can be given 
as 

A/15 = ^TrE {{pE Ps) HintHint " 2Hint {PE ® Ps) Hint + HintHint {PE <^ Ps)} 

- \ {Ps{'Hint){nint) - 2(Hint)/is(Hint) + (Hint) (Hint) Ms} (B9) 
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An interaction Hamiltonian can always be written as 

^max 

"^int ^^fi<» 

i 

where fi are Hermitian operators on the environment Hilbert space, while ipi are Hermitian operators on the system 
Hilbert space. With this form of the interaction in mind, one can write compactly 



A/is = 



I E - (/')(/])) {^Mt^s - 2^]^^s^, + fls'P^v]} 



-^E((/4)-(/^)(/l)){ 



] 



(BIO) 



where, as before, (...) = Tr£;(p£; . . . ). ity of f's. As a result, 

^From its definition jBTj, operator A satisfies 2 A = 

« ('^)o + (f )o_il^«^ look at the matrix elements of dC,s „ , I ( dg\ ^^^^^ 

the matrix [Eq. ^] -dT^^' " 2 I i ' ^^^^^ 



{hf])-{h){f]. 



Furthermore, m is a positive-definite matrix, since for 
we immediately note that = m*^ due to the hermitic- any complex imax column c— (ci, . . . , Ci^„^ it holds 



{c\m\c) = 5]c"™.,c^" = ^c" {(/,/]) - {h){f])]d 



(B12) 



Tr^ < PE 



Trs < /9_E 



> 0. 



The last inequality can be obtained by noting that, for though, unlike the Lindblad dissipator that has the units 



any matrix a, 

Tr^ [pEaa)) 



k,k' 



> 



Ek I k 
Pe Ffc 



> Y.^P%f\al\'> 

k 

= iTrfi (psa)!^ . 



Y.PE^'k 



(B13) 



As a result, we conclude that —A has the form expected 
from the Lindblad dissipator (it has the units of t~^. 



J 



C12M can be found as 



of <-i). 

Up to the second order in time, the generator /C of the 
non-Markovian map (jl3[) can now be approximated as 

/C(t) « -iCaS - 2At + o{t^), (B14) 
where CcS = £off(0) from Eq. (jB5ap . and A is given in 
Eq. (lB8|) . 



1. How to calculate £i2A4 and A 



dE-l 



a'P'\a'f3' 



(B15) 
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For a fixed a, (3, there is a d^-dimensional space spanned by all \ia, if3). A unit operator in this space can be written 
3 J2i=i^ \bi,ai3) {bi,ai3\ + \af3) {a(3\ = lap 



, 

V (^\L\b,,a'p') - (^|L {la'p' - (^1) \a'(3') 



a/3|L|a'/3'; / 



dE 



i— 1 ^ 

dE 

[(^int)iQ''5;^/ - (Hint)-;^ ^2' 



The last line is easily obtained by showing that the contributions from the environment Hamiltonian [(Ticnv)^^ = 
(Hb)-;^^] and from the system Hamiltonian [("Hgys)-^ = (Hs)^] vanish. 

When one deals with interaction Hamiltonians of the hopping type, i.e., those that contain an odd number of 
environmental creation/annihilation operators and therefore necessarily alter the environmental state, all (Tiint)^"' — 
0, and clearly LnM. — 0, which we used in Sec. [C] Also, when the statistical operator is uniform (p^ = 'Pe)i 
L\iM. — 0. Note how this term accounts for the information influx from the environment, because it captures the 
deviation of the environment statistical operator from the uniform statistical operator (the uniform statistical operator 
carries the maximum information entropy, i.e., environment has no information to transmit). 

In order to calculate A, which was defined as 2A ~ C12C21 + C.\2C.i2M. — Lx^M. (£11 + £12^!) in Eq. (jB7p . we 
should first note that LyiM. (£11 + L\2M) is commutator generated, i.e., 

UriM (£11 + lix2M) - [(Hint) - Hint, {Us + (Hint), ..•]]■ 

The term £12/121 + £i2£22-A^ can be rewritten as 

(£i2£2i+£i2£22X)::f^, = (^\L' ~ LPL\^') - (B17) 

{a'(3'\a'f3') 

ds 



(aPlL^la'/]') ^ - V {a(3\Lm){ja\L\a'(3')^^^=^ 



Y^{ia,i(3\L^\]a',3(5')fy'E 

^ dm ds 

— ^ ^ p^(ia,z/?|L|/c7,fccr)(fc7,A:cr|L|ja', j/3'), 



i,j,k=l 7,o'=l 

where the eigenbasis of the environment initial statistical operator pE is chosen to be the environmental basis. Upon 
a straightforward (and somewhat lengthy) calculation, with the only constraint being that [pet'He] = 0, which is 
typically satisfied, we obtain Eq. (|B8p . 

I 



APPENDIX C: TWO ADDITIONAL EXAMPLES 

The following two examples serve to illustrate that the 
usefulness of the coarse-grained map (j22p may extend be- 
yond the strict validity specified by p4|) . and may offer a 
particularly simple way to identify the steady state alone 
from first principles. 

The first example (IC 1|) is analytically solvable and pos- 



sesses the long-time Markovian evolution regardless of 
the interaction strength. We show here that there exists 
a mathematical coarse-graining time r, shorter than any 
other timescale in the system or environment, so that 
the exact long-time Markovian evolution coincides with 
that obtained from the short-time evolution by coarsen- 
ing over T (|22p . 

On the second example (|C 2p . we show that relaxation 
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towards the correct equilibrium state is easily obtained 
by using (|22p (or equivalcntly by employing Theorem 2 
in Sec. 1111 . 



1. Spin-boson model with pure dephasing 

One of the few analytically solvabl e ^^'^"^i^^'^^i^^'^^ open 
system problems is that of a two-level system coupled to 
a dephasing-only boson bath, with the relevant Hamilto- 
nians are given by 

Us = ^a.,HE^Y. ^1 ^) ' 

Hint = ^a.{<?(r!,-)6,- + 5(%)*&^4- ^^^^ 



Here, cr^ is the Pauli matrix, bt and feg- and the boson 
creation and annihilation operators of the q-th boson 
mode, respectively, ±co'/2 are the system energy levels 
(divided by h), and fi^ is the boson mode frequency. 
The boson modes are initially in a thermal state with 
(ug) = {h\bq) = exp(;if2,/fc^r)-i - Because of the inter- 
action linear in environment creation/annihilation oper- 
ators, (TYint) ^ 0, SO Cs = CcS'- 




10 
0-10 




(C2) 



where the rows/columns are ordered as 1 = |+) (+| , 2 = 
|-K)(-|,3 = = |-)(-| (± refer to the posi- 

tive/negative (upper/lower) energy state). Operator A 
can be calculated according to (jBSP as 



A 





10 

10 





2i:wn,-)P».i.(^) 

q 

2 dl]p(f7)|g(f7)pcoth(^^j. 



(C3) 



where I?(ri) is the density of boson states. 

Cs and A obviously commute, and their common zero 
eigenspace [A/'(A) — A/'(£off)] contains all density matri- 
ces with zero off-diagonal elements. This means that, for 
a given initial statistical operator, the off-diagonal matrix 
elements will decay to zero while the diagonal elements 
remain unchanged: 



[Ps{t)]++ 

[P5(i)] + _ 



[PS(0)]++, [p5(0]-_ 

[p5(0)]_+e+'"*-^''^*, 



(C4) 



The steady state will be determined by simply annulling 
the off-diagonal elements. This is the correct steady 
state, as shown in the exact solution^^. 

Instead of exp(— A^Tt), in the exact solution decoher- 
ence is seen through the term exp[— r(t)], where the r(i), 
the decoherence exponent, behaves as 



m = f 

Jo 



dn2V{n)\g{n)f coth 



/ nn \ sin2(m/2) 



{n/2y ■ 

(C5) 

For short-times, T{t) « A^t^, as should be expected, 
because we know our expansion is exact up to 

the second order in time. In the long-time limit for 
r(t), only the low frequency contributions survive, since 



lim. 



sin"'(Ot/2) _ 
(0/2)^t 



SO 



oo 



t hni 27rX'(r2)|5(r2)pcoth 



2kBT 



(C6) 



We need to match this long-time behavior of r{t) with 
our coarse-grained term XdTt, in order to obtain r. 



lim 

t^oo Arfi 



limo^o 2nV{n)\g{Q)\^ coth (^1^) 
2/o°°df7P(f7)|g(f7)Pcoth(5^) 



(C7) 

Let us consider the example of an Ohmic bath (e.g., page 
228of Ref. 38,), with D{n)\g{n)\^ = il7 exp(-0/17c) and 
flc being a density-of-states cutoff frequency. Typically, 
Mlc 3> fc^T. In the numerator, one can approximate 

(2!^) ~ ^Sr^ while the coth function in the de- 
nominator is always greater than 1, yielding 



/kBT\ 27r 2n 



(C8) 



Being typically the largest frequency scale in the full SE 
problem, fie sets the shortest physical timescale. Clearly, 
T is even shorter than the period associated with Qc, 
which justifies our use of the short-time expansion and 
subsequent coarse-graining. 

Note the long-time behavior exp(— t/r^) of the deco- 
herence term F, where tt — h/nkBT is the thermal cor- 
relation time. However, our time r is the mathematical 
coarse-graining time, which is very short. The relation- 
ship between the correct physical correlation loss time 
and the mathematically appropriate time is 



(AdTr) 



(C9) 



2. Jaynes-Cummings model in the rotating wave 
approximation 

The Jaynes-Cummings Hamiltonian in the rotating- 
wave approximation22iii22iii2iii2^ describes the decay of a 
two-level system in the presence of a single boson mode 
of resonant frequency. The relevant Hamiltonians are 
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1 



Hint = g {b^a- + ba+) 



(ClO) 



Here, Cz, cr+ — [ox + iCy) /e, and ct- = ((To; — i(Ty) /2 are 
the Pauli matrices, 6^ and 6 are the boson creation and 
annihilation operators, respectively, ±ti;/2 are the system 
energy levels (in units of frequency) and uj is also the 
boson mode frequency, and g is a parameter measuring 
the interaction strength. The boson mode is initially in 
a thermal state with (n) — {Wb) — c:K.p{huj}kBT)-i ■ 
in the spin-boson example, L^s = because of the 
interaction linear in environment creation/annihilation 
operators: 



-off 




10 
0-10 




(Cll) 



Operator A can be calculated according to Equation 
B8t as 



2(n)+ 2 ~2{n) 

2{n) + 10 

2{n) + 1 

-2{n) - 2 2{n) 



(.C12) 



Cs and A commute, and we immediately note two com- 
mon one-dimensional eigenspaces: p-\ is associated with 



the Cs and A eigenvalues u and g'^{2{n) + l)/2, respec- 
tively, while p (- is associated with the eigenvalues —uj 

and g'^{2{n) + l)/2. 

On the other hand, the space spanned by and 
|— )(— I is the null space of £off- Solving the eigenproblem 
of A reduced to this space gives 



det 



gmn) + l)-X -g^n) 
-gH{n) + l) g'{n)-X 

A = and A = = c/^ (2(n) + 1) . 



(C13) 



An eigenvector pP = (/^5J__|_, corresponding to the 

zero eigenvalue of the matrix A is characterized by 



1 



(C14) 



If we are looking for a statistical operator that belongs 
to the zero cigcnspacc of A, it also has to satisfy the 
constraint of the unit trace, which fixes 



2{n) 



1 



(n) + 1 



2(n) 



1 



(C15) 



One recognizes these components as the thermal equilib- 
rium values of the population of the upper and lower level 
of our two-level system, respectively (see, for instance, p. 
149 of Ref. Issl ). Therefore, by seeking the steady state in 
A/'(A) n A/'(£cff), we have obtained the physically correct 
result. 
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